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ABSTRACT 

We use the science demonstration field data of the Herschel- ATLAS to study how 
star formation, traced by the far-infrared Herschel data, is related to both the 
accretion luminosity and redshift of quasars selected from the Sloan Digital Sky 
Survey and the 2SLAQ survey. By developing a maximum likelihood estimator 
to investigate the presence of correlations between the far-infrared and optical 
luminosities we find evidence that the star-formation in quasar hosts is correlated 
' with both redshift and quasar accretion luminosity. Assuming a relationship 

of the form L m oc L e QSO (l + z) c , we hnd 9 = 0.22 ± 0.08 and C = 1-6 ± 
0.4, although there is substantial additional uncertainty in £ of order ±1, due 
to uncertainties in the host galaxy dust temperature. We find evidence for a 
large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic 
dispersion in the correlation between Lqso aud £ir, suggesting that the latter 
may be due to a direct physical connection between star formation and black 
hole accretion. This is consistent with the idea that both the quasar activity 
and star formation are dependent on the same reservoir of cold gas, so that they 
are both affected by the influx of cold gas during mergers or heating of gas via 
feedback processes. 
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1 INTRODUCTION 



AGN activity is now widely believed to be an important 
phase in the evolution of every massive galaxy in the 
Universe. This belief stems firstly from the discovery of a 
relatively tight correlation between the mass of a galaxy's 
bulge and th e mass of its central supermassive black hole 
(SMBH: e.g. iMagorrian et all 1 19981 : iFerrarese fc MerrittJ 
l200d ; iGebhardt et alj l200Ch which implies that build- 
up of the stellar mass and t he SMBH mass are cau sally 
connected (although see also lJahnke fc Maccioll20ld . who 
show that averaging of parameters by mergers can also 
produce such a correlation). 

From a more theoretical perspective, AGN have 
become a key ingredie nt for semi-analytic models of 
galaxy formation (e.g. Benson et alj 20031; Cole et all 



l200d : iGuiderdoni et al J 1 19981 : IGranato et al.ll2o"oil ). The 

over-production of stars at the bright end of the galaxy 
luminosity function by these models has led to the inclu- 
sion _of_aii_extrj^mirce of heating or "AGN feedback" (see 



e.g. iBower et al]|200d ; ICroton et aTll200d ; IGranato et aH 
2004). However, it is still unclear what kind of AGN- 
driven feedback is the most important, e. g. two com- 
pleme ntary mechanisms are described by ICroton et al.l 
(2006): "quasar-mode" feedback, proposed to be effec- 
tive during episodes of efficient accretion of cold gas, and 
"radio-mode" feedback, caused by the relatively ineffi- 
cient accretion of hot gas (in the absence of cold gas) 
when an AGN would have lower optical luminosity, but 
would also be generating radio jets. Thus, it is now clear 
that if we are to understand galaxy formation and evo- 
lution, we must also obtain a clear picture of AGN, their 
environments and the co-evolution of SMBHs and their 
host galaxies. 

It is difficult to study the host galaxies of dis- 
tant quasars, due to the nuclear quasar emission over- 
whelming the stellar emission at optical and near- 
infrared wavelengths. However, Hubble Space Telescope 
and high spatial-resolution ground-based imaging have 
enabled the inve stigation of the stellar populations of 
quasar hosts (e.g. Bahcall et al . 1994; McLure et al .1 19991 ; 
iDunlop et alj|2003l ; iKotilainen et alj 12009 ). demonstrat- 
ing that most luminous quasars reside in massive el- 
liptical galaxies. A number of studies have attempted 
to determine the star-formation properties of quasar 
host galaxies, with some using optical colours (e.g. 



Sanchez et alj 2004h or spectroscopy (e.g . Nolan et alj 



200 ll ; ISilverman et all [20091 ; iTrichas et at] l201Ch . Some 
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information on the host galaxies of AGN has been 
gleaned from powerful radio galaxies, where the quasar 
nucleus is obscu red due to the putative dusty torus (e.g. 
lAntonuccll Il993l ). Again these studies sho w that pow- 
erful AGN reside in massive galaxies (e.g. Jarvis et al] 
l200ll ; iMcLure et ai]|2004l ; ISevmour et alj|2007h and there 
appears to be a trend where the most powerful exhibit 
signs of recent star formation fe.g. iBaldi fc Capettill2008l ; 
iHerbert et al]|20ld ). However, optical colours and emis- 
sion lines can be contaminated by the quasar, so that 
its contribution must be modelled to estimate the star- 
formation rate, and all of these studies suffer from the 
uncertainty arising from having very little information 
about how much star formation is obscured by dust, and 
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thus how much of the emission due to young stars is re- 
processed to far-infrared and sub-millimetre wavelengths. 



Spectroscopy with the Spitzer Space Telescope has 
enabled the determination of star-formation rates in in- 
dividual AGN via mid-infrared emiss ion lines, which 



are much less susceptible to dust, (e.g. Schweitzer et al.| 


20061: Netzer et al. 


120071: Martnicz-Sansigre et al.l 


2008|; 


Lutz et al. 


2008; 


Shi et all 20091: Trichas et all 


2009: 


Lutz et al. 


120101). However, the relatively low sensitivity 



of Spitzer's infrared spectrograph (IRS) has meant that 
only small, severely flux-limited samples of objects could 
be investigated, making the precise relationship between 
star-formation rate and quasar luminosity difficult to ro- 
bustly determine. 

(Sub) millimetre instruments such as SCUBA and 
M AMBO have provided important information on repro- 
cessed emission from star formation in distant quasars 
and radio galaxies, e.g. photometric surveys revealing lu- 
minous (~ 10 13 L Q ) far-infrared emission from a number 
of radio galaxies and quasar hosts, making it clear that 
major (SFR ~ IOOOMq yr _1 ) episodes of star formation 
can c oexist with powerful AGN (e.g. [Archibald et al. 
200l |: IPriddev et al.ll2003l:IOmont et alJliool : IPaee et all 
20041 : iMartmez-Sansigre et al.ll2009l v" Furthermore, sub- 
mm imaging has revealed overdensities of sources in the 
vicinity of individual high-z AGN, plausibly belonging 
to the same large-s c ale structures llStevens et al.l 120031 : 
IPriddev et all 120071 : IStevens et all l20ld ). Conversely. 
cross-correlating SCUB A data with deep X-ray surveys, 
I Alexander et al.l |2005) also show a large fraction of 
sub-m m galaxies contain buried AGN, while lLaird et all 
(|201(]| ) find that the X-ray AGN fraction in sub-mm 
galaxies is 20-29 %. 

The main drawback of the sub-mm studies to date 
is that, due to the high background associated with 
ground-based observing and the small number of pix- 
els in existing instruments, it has not been possible to 
obtain very large samples of obje cts. A recent study by 
ISerieant fc Hatziminaogloul |2009l . hereafter SH09) at- 
tempted to overcome this shortfall by combining data 
from a wealth of Spitzer and ISO observations at mid- 
and far-infrared wavelengths to investigate how obscured 
star formation in the hosts of quasars evolved over cos- 
mic time and also how this may be related to the bolo- 
metric output of the AGN. They found evidence for an 
increase in far-infrared luminosity as a function of both 
cosmic epoch and quasar luminosity. However, they mea- 
sured the correlation between far-infrared and quasar lu- 
minosity using wide redshift bins, which can potentially 
confuse the effect of redshift with that of quasar luminos- 
ity, since the latter evolves strongly with redshift in flux- 
density-limited samples. In addition, their work extrapo- 
lated the far-infrared flux from mid-infrared wavelengths 
using SED templates, and it would clearly be preferable 
to measure the far-infrared flux directly. 

The recently-la unched ESA Herschel space observa- 
tory fPilbratt 2010) offers an unprecedented combination 
of sensitive detectors, large (3.5-metre diameter) collect- 
ing area, and low background, which makes it possible, 
for the first time, to conduct deep surveys of large areas 
of the sky at far-infrared/sub-mm wavelengths. 

In this paper we present a study of the far-infrared 



fluxes of known quasars in science demonstration obser- 
vations of the Herschel Astrophysical Ter ahertz Large 
Area Survey (H-ATLAS: lEales et all |2010| V This is the 
first single sample of quasars measured in the sub-mm 
which is large and diverse enough to directly determine a 
relationship between far-infrared (cool dust) luminosity 
and optical (quasar) luminosity, without using infrared 
information at wavelengths short enough to be possibly 
contaminated by wa rm dust in the AGN torus. While 
ISerieant et al.l (|201Cf ) determine a relationship between 
optical and far-infrared luminosity using this data in com- 
bination with infrared luminosities derived from shorter- 
wavelength Spitzer and ISO data, here we use only the 
H-ATLAS data. 

Section[2]describes the H-ATLAS observations used. 
Section explains how we construct our quasar sam- 
ple, determine far-infrared fluxes for them, and find a 
maximum-likelihood fit to the relationship between far- 
infrared luminosity, optical luminosity, and redshift. The 
results of this fitting are discussed in Section U and con- 
clusions are presented in SectiorO 



2 OBSERVATIONS 

When completed, H-ATLAS will provide photometry for 
550 square degrees of the sky in five far-infrared/sub- 
mm bands. The 5-er depths at wavelengths of 100, 160, 
250, 350, and 500 fj,m are 132, 126, 32, 36 and 45 mjy 
respectively. 

In this paper, we use only data from HerscheVs 
SPIRE instrument (jGriffin et all 120101 ). at 250, 350, 
and 500 fj,m, s i nce t he PACS 100 and 160 |im data 
(|Poglitsch et al.l |2010T ) are not deep enough to provide 
useful constraints on the infrared lumino sities of the 
quasars we are studying (|lbar et aljkoich . The longer 
SPIRE wavelengths are also less likely to be contami- 
nated by direct emission from the warm dust of the AGN 
torus. This means that the emission can be more reli- 
ably attributed to cool dust heated by star-formation 
processes, and its spectral energy distribution (SED) ap- 
proximated by a simple greybody. 

This paper is also restricted in area to the H-ATLAS 
science demonstration (SD) field, which covers approxi- 
mately 16 square degrees - around 3 per cent of the total 
survey area. This SD data has been publicly released, and 
both the maps and catalogues can be retrieved from the 
H-ATLAS webpageQ. 



3 METHOD 

Our quasar sample consists of all objects in the field 
which are identified as quasars i n the seventh data re - 
lease of the SDSS Quasar Catalog (|Schneider et al.ll2010l L 
or the 2dF-SDSS LRG (luminous red galaxy) and Qu asar 
spectroscopic catalogue (2SLAQ;|Croom et al.ll2009h . Af- 
ter removing the known blazar HB0906+015, since its 
far-infrared emi ssion will not be dominated by star for- 
mation (see e.g. lGonzalez-Nuevo et ai1l2010l ). we find 372 

1 http://www.h-atlas.org 
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quasars in the whole SV field, with a reasonably conser- 
vative cut on the area (by eye) to avoid using the noisier 
data at the edges of the field. 

Of the 372 quasars, only 29 ha ve matches to ob- 
jects in the 5-0- H - ATLAS catalogue (|Rigbv et al.ll2010l : 
ISmith et alj|20ld ). To improve our statistics, we mea- 
sure fluxes di r ectly from the SPIRE maps described by 
IPascale et all (j2010h at the position of each quasar. To 
do this, we take a point measurement of the flux in the 
nearest pixel to the SDSS position on each PSF-filtered 
and background-subtracted map. 

To account for source confusion, the most likely flux 
from the object (which is used to produce Figure[2]) is ob- 
tained by subtracting the mean pixel value in the map. 
However, when we perform a maximum-likelihood fit for 
the coefficients of Equation [5] as described below, there 
is an implicit assumption that the measurement noise is 
symmetric about zero. The confusion noise is very sig- 
nificantly skewed, and can cause false correlations to be 
found in random data, so for the fluxes used in these fits 
we symmetrise the noise distribution by subtracting the 
flux in a different randomly selected pixel from each flux 
measurement. This increases the overall noise level, so we 
average the likelihoods from a number of random back- 
ground realisations in order to reduce the noise level in 
the final result. 

We also construct a comparison sample of fluxes 
drawn from random positions in the maps, with their 
backgrounds subtracted in precisely the same way as 
the real fluxes, and are then used to produce maximum- 
likelihood estimates in the same way as the real data. 
Comparing these estimates with those from the real data 
allows us to increase our confidence that any trend de- 
tected in the data is due to the influence of the quasars, 
and not merely a result of the noise properties of the 
maps, which have a significant impact due to our inclu- 
sion of individual noisy objects instead of using binned 
data. 

To compare the AGN luminosity of the quasars with 
that from star formation, we make the assumption that 
the optical light is dominated by the AGN, while the far- 
infrared/sub- mm is dominated by the cool dust heated 
by star formation. This latter assumption i s con sistent 
with the results of lHatziminaoglou et al.l ([2010) , who 
fit models with AGN and star-formation contributions 
to the mid- and far-infrared emission from AGN hosts 
and find that star formation dominates the far-infrared 
fluxes. Their result depends on the assumption that the 
dusty torus heated directly by an AGN has a fairly sim- 
ple, approximate ly toroidal geometry, as described by 
iFritz et alj (|2006t ). The infrared SEDs of AGN can, how- 
ever, be full y reproduced by mor e complex clumpy torus 
models (e.g. iNenkova et"aL l2008t ). so, at least until these 
models are tested by high-spatial-resolution imaging of 
AGN tori, we cannot completely rule out the possibil- 
ity of a contribution to far-infrared fluxes from dust 
heated directly by the AG N. However, the fact that 
lHatziminaoglou et alj (|2010h find that the SPIRE colours 
of AGN are indistinguishable from those of normal galax- 
ies suggests that the AGN does not make a strong direct 
contribution at these wavelengths. 

All of the quasars have spectroscopic redshifts from 
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Figure 1. Absolute i-band magnitudes (calculated as de- 
scribed in the text) for our sample of spectroscopically con- 
firmed quasars, plotted as a function of redshift. Objects with 
spectra from SDSS are shown as black diamonds; those with 
2SLAQ spectra have red squares. 



their respective surveys, and optical photomet ry in the 
ugriz bands from the SDSS (|York et alj|2000h . We use 
the se data to estim ate the quasar luminosity, follow- 
ing ISchneider et al.l and others in using the absolute i- 
band magnitude, Mi, as a proxy for the luminosity of 
the central engine. This is calculated from the observed 
i magnitude (after the subtraction of Galactic extinc- 
tion) and the redshift, by assuming a power-law SED 
(S v oc v a ) with spectral index a — —0.5 and assum- 
ing a ACDM cosmology with Q m = 0.3, Qa = 0.7, and 
Ho = 70kms _1 Mpc" 1 . Figure [T] shows the distribution 
of optical luminosities of our sample as a function of 
redshift, indicating the source of the spectroscopic red- 
shift (SDSS or 2SLAQ, or both). From this we can see 
that the 2SLAQ measurements in this field allow the 
inclusion of quasars at least one magnitude fainter (at 
z < 3) than the SDSS spectroscopic survey probes. This 
greater range in optical luminosity is crucial for decou- 
pling optical luminosity and redshi ft, and gives our study 
an ad vantage over similar work by lHatziminaoglou et al.l 
(2010). We note, however, that the range of bolometric 
luminosities spanned by our quasar sample is 10 45 erg 
s _1 < I/boi < 10 48 erg s _1 , so all of the objects are firmly 
in the quasar (as opposed to Seyfert) regime. Thus, if 
there are two "modes" of accretion with different feed- 
back properties, this study will only be relevant to un- 
derstanding the "quasar mode" , making it u sefully com- 
plementary to the work of IShao et al.l |2010l ). which uses 
different Herschel data to probe the star-formation rates 
in lower-luminosity X-ray-selected AGN. 

As mentioned earlier, the SPIRE photometry at 250, 
350, and 500 /im is unlikely to be contaminated by the 
AGN directly for even the higher redshift quasars in the 
sample, and is deeper than the H- ATLAS PACS photom- 
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Figure 2. Estimated total infrared luminosities (calculated 
by integrating the best-fitting greybody with T = 30 K and 
P = 1 .5 over the 8-1000 fim wavelength range) for our sample 
of spectroscopically confirmed quasars, plotted as a function of 
redshift. Objects with spectra from SDSS are shown as black 
diamonds; those with 2SLAQ spectra have red squares. Error 
bars are plotted where Lir is greater than the uncertainty 
in the estimate, otherwise 1-<t upper limits are shown using 
arrows. 



etry jlbar et al.ll20ld ). so we use only the SPIRE fluxes to 
constrain the overall infrared luminosity. To make the re- 
sults directly comparable with each other, we fit a simple 
isothermal grey dust model, with flux density 



S„ oc ^/(e^W - 1) 



(1) 



where h and k are Planck's and Boltzmann's constants, 
v is rest-frame frequency, and where we assume fixed val- 
ues for both temperature T — 30 K and ft = 1.5. We then 
integrate the greybody curve between rest-frame wave- 
lengths of 8 and 1000 /mi to obtain an estimate of Lir. 
(This wavelength range is chosen for convenience and for 
comparison with other studies, and is not intended to 
imply that the cold dust represented by the greybody 
dominates the SED at wavelengths as short as 8 /mi.) 
We estimate the uncertainty on Lir by increasing it un- 
til the x 2 value for the fit to the data increases by unity. 

Since we include far-infrared flux measurements at 
our quasar positions regardless of whether there is a sig- 
nificant detection in the far-infrared, fitting the normali- 
sation of a fixed template is the only robust measurement 
possible. This procedure does not account for the possi- 
bility of dust temperature or emissivity evolution either 
with redshift or quasar luminosity, which may well be im- 
portant effects to consider. However, we are able to use 
the well-detected objects in our sample to address the 
suitability of our choice of T and ft. 

To check our choice of greybody temperature, we fit 
greybody SEDs to the SPIRE fluxes of individual quasars 
with good detections (signal-to-noise > 3) in all three 



SPIRE bands, with ft fixed at 1.5 and T allowed to vary. 
The best-fitting temperatures for eight out of the nine 
well-detected quasars were very similar, with mean 32.4 
K and a standard deviation of 2.8 K, while one outlier 
fits a temperature of 58.5ly 4 g° K. The eight quasars span 
a redshift range of 1.1 < z < 2.3 and a range in infrared 
luminosity (estimated by integrating the greybody be- 
tween 8 and 1000 /mi) of 6 x 10 12 < L m /L Q < 2 x 10 13 . 
They are consistent with the temperatures of g alaxies as 
a whole as determined bv lAmblard et all l|2010l l. who find 
a mean temperature of 28 ± 8 K for galaxies detected at 
a significance of 3-<r in the H-ATLAS SDP field, and a 
mean temperature of 30 ± 9 K for a compilation of galax- 
ies in H-ATLAS and a number of other far-i nfrared sur- 
veys. This consistency is unsurprising, since lElbaz et all 
(|20ld ) also found (using deeper data from Herschel) that 
the dust temperatures of radio-quiet AGN hosts were not 
significantly different from the normal galaxy population. 
Thus we consider that a temperature T = 30 K is likely 
to be a reasonable choice for our quasar sample, but is 
not completely certain, and we return to this point later. 

The outlying object has substantial uncertainties 
on its temperature, so we do not consider it an indi- 
cation of a serious problem with our use of a single 
temperature for all objects. However, the fact that it 
is at a higher redshift (z = 3.4) and slightly higher 
luminosity (Lir/Lq = 3 x 10 13 ) than the other ob- 
jects means that it is broadly consistent with the evo- 
lution of dust temperature, to higher temperatures at 
higher redshif ts and luminosities, s een in several pre- 
vious studies ([C hapman et al.l 1200 5; Kovac s et al. 2006; 
iMichalowski et al.ll2010l : lAmblard et alj201Ch and others . 

We do not attempt to use a luminosity- or redshift- 
dependent dust temperature in this work. Increasing the 
dust temperature serves to shift the greybody peak to 
shorter wavelengths, which would result in a higher lu- 
minosity fit to our data (which are mostly longward of 
the peak) . Thus, if the temperature of the dust in quasar 
hosts does increase with redshift and/or luminosity, then 
we will systematically underestimate Lir for objects at 
high redshift and Lir. However, since the trend in mean 
temperature found by lAmblard et alj is weak compared 
with the dispersion around the relation, we do not expect 
this to be a very significant effect. 

To get an indication of the effect of temperature (and 
its uncertainty) on the parameters we derive, we perform 
fits to all objects using both T = 30 K and T = 45 K, 
where we have selected 45 K because it has a similar peak 
wa velength to the averag e z > 4 quasar SED constructed 
bv lPriddev fc McMahonl (|200lh (who find parameters of 
T = 41 ± 5 with ft = 1.95 ± 0.3). We also note that, 
while ft may not in fact be 1.5, the luminosity estimates 
obtained using ft = 2 (using the best-fitting temperature 
at this ft) differ by less than 1 per cent, so we consider 
that the selection of a value of ft is broadly degenerate 
with the choice of T, at least for the purpose of estimating 
infrared luminosity from SPIRE fluxes. 

Following SH09, we assume that the dependence of 
the infrared luminosity on the optical quasar luminosity 
and redshift can be expressed as power laws, i.e.: 



Li 



al: 



QSO 



(1 + *) C 



(2) 
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with intrinsic (Gaussian) dispersion in 8 and C of A8 and 
AC respecti vely. 

Unlike ISerieant fc Hatziminaogloul . who stack in- 
frared data to look for trends in bins of redshift and 
optical luminosity, we fit for the power law coefficients 
of equation p]l. and the normalisation A, by finding the 
peak in the likelihood distribution calculated over all data 
simultaneously. Our approach has the substantial disad- 
vantage that variations in the power laws as a function 
of redshift or luminosity will be hidden from us. How- 
ever, under the assumption that the values of C and 8 are 
the same at all redshifts and luminosities, our method 
enables us to make the most sensitive possible determi- 
nation of these parameters, since we do not lose informa- 
tion through binning. Additionally, we are immune to the 
potentially confusing effects of correlations within bins. 

We incorporate the dispersions A8 and AC by com- 
bining them with the uncertainty <Ji on each data- 
point. The likelihood, Prob(data|0, C, A8, AC, ^4) at a 
given point in parameter space is then given by 



Prob(data|0, C, AO, AC, A) = T] 



where 

2 _ (ilR,data,i — ^IR, model, i ) 2 
Xi 



"X?/2 



V^K 2 + ° l 



e,i 



+ cr, 



:(3) 



o~e,i = L m , model, i In Lqscm A6> , and 



= Lt 



,iln(l + z) AC 



(4) 

(5) 
(6) 



The prior space on the parameters 8, C, A8, AC, and 
A, over which the likelihood was evaluated, was chosen 
iteratively to enclose more than 99.999 per cent of the 
likelihood. The priors were linear in 8, C, A8, and AC, 
and logarithmic in A, with the limits A8 and AC ^ 0. 



4 RESULTS AND DISCUSSION 

Figure [3] shows the results of the likelihood evaluations, 
marginalised over all parameters except the two shown in 
each panel. This is the result of averaging the likelihoods 
obtained in 17 random realisations of the background, as 
described in Section [2] Contours labelled 1-a, 1-a, 3-ct, 
and 4-cr enclose 68.3, 95.4, 99.73, and 99.994 per cent of 
the total likelihood. We show the results of fits assuming 
both T — 30 K and T = 45 K to give an indication of 
the way in which our uncertainty in T should contribute 
to the uncertainty in the derived parameters. 

This figure shows that, assuming the dust in quasars 
has a temperature T — 30 K, there is good evidence for 
correlations of Lir with both Lqso and z, since the max- 
imum likelihood values of both 8 and C are significantly 
non-zero, and that while higher temperatures can mimic 
the effect of the correlation with z, the correlation with 
Lqso is robust to simple changes in the overall dust tem- 
perature. For the T — 30 K case, the best-fitting values of 
each parameter, with their 68.3 per cent confidence lim- 
its (marginalising over all other parameters in each case) 
are: 



c 

AO 
AC 



= 0.22 ± 0.08 

= 1.6 ± 0.4 

= 0.08 ± g;£ 

= 3.4 ±g : ? 



(7) 
(8) 
(9) 
(10) 



Although, as can be seen in the lower left panel of 
Figure [3J there is some degeneracy between 8 and C, our 
data exclude 8 = at approximately the 2.5-cr level. The 
(flat) prior on 8 allowed it to take values between -1.5 
and 1 (i.e. a much larger range than we show in Figure 
O, so we are confident that this is not simply an artefact 
of reaching the edge of the prior range. 

The other panels of Figure [3] show that there is not 
a strong degeneracy between any other parameters, so 
it is not possible to reproduce the effect of a correlation 
of LiR with either Lqso or 2 by simply increasing the 
intrinsic dispersion in the other parameter; indeed these 
panels suggest that approximately the same result would 
be obtained if dispersion in the power law indices was 
ignored altogether. 

While the figure shows that the detected evolution 
of Lir with redshift is very strongly dependent on the 
assumed value of T (although a temperature as high as 
45 K for the whole sample seems unlikely) , the correlation 
with Lqso seems to be robust to this uncertainty. We find 
a slightly weak er inde x, 8, for the correlation with Lqso 
than found by SH09, whose average result (8 = 0.44 ± 
0.07) is shown in figure [3] as the point with errorbars. 
This may in part be due to the relatively large redshift 
bins (e.g. 0.5 < z < 1, 1 < z < 2) which were used in 
SH09, since, because of the strong degeneracy between 
8 and C, these could allow trends with redshift to affect 
objects within a single bin. However, the disagreement is 
not very significant. 

In contrast, performing the same maximum- 
likelihood fit to data with far-infrared fluxes drawn from 
random positions in the maps (with the same position 
used at each wavelength, and the dust-fitting performed 
assuming a temperature of 30 K) we obtain the results 
shown in Figure [4] The important thing to note is that 
the best-fitting value of 8 is zero for the random flux case, 
indicating that the result we have obtained is not simply 
due to the noise properties of the far-infrared maps, but 
rather is due to far-infrared flux genuinely associated with 
the quasars. (Note that the non-zero value for C derived 
from the fit to random data is not unexpected, since the 
real quasar luminosities and redshifts were still used in 
this case, and these are themselves correlated.) 

Using deeper far-infrared data from Herschel, but 
with a narrower range in optical luminosity (because 
only SPS S, and not 2SLAQ, quasars are available in 
the field) , lHatziminaoglou et all (|2010h find a very sim- 
ilar result, 8 = 0.35 , for objects at redshifts z > 2. 
ISerieant et all (|2010T ) use a compilation of H-ATLAS 
and other data to find a weakly declining correlation 
between Lir and quasar luminosity, parametrised as 
8 = (0.5875 ± 0.045) - (0.09275 ± 0.012) z for redshifts 
< z < 4, which is a slightly higher value of 8 than we 
find, but broadly consistent within the errors and at the 
redshifts we are sensitive to. Our value of 8 is also consis- 
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Figure 3. Likelihood contours for pairs of the parameters 8 and f (from equation [2} and their intrinsic dispersions A9 and A£, 
marginalised over the parameters not displayed in each panel. Contours labelled 1-<t, 2-cr, 3-cr, and 4-cr enclose 68.3, 95.4, 99.73, 
and 99.994 per cent of the total likelihood. The poi nt wit h error bars in the lower left-hand panel shows the average values of 
9 and f, and their 1-tr uncertainties, determined by SH09 via stacking in the mid-infrared. Solid contours show a fit to objects 
assuming T = 30 K; dotted contours show a fit to the objects assuming T = 45 K. 



tent with the work by ISilverman et "all (2009), who find 
9 — 0.28 ± 0.22 for the slope of the power law relating 
[On] luminosity (tracing star formation after correcting 
for AGN contamination) to X-ray luminosity (tracing the 
AGN). Silverman et al. study a sample of AGN over a 
narrower range of redshifts (0.5 < z < 1) and with some- 
what lower luminosities than most of our sample (10 44 ' 3 
erg s _1 <Lboi~ 10 45 ' 8 erg s _1 , assuming a bolometric cor- 
rection I/boi/£2-iokcV ~ 30); the consistency of their 
result with ours suggests that the relationship between 
Lib, and Lqso ma y not vary with LQSO as suggested by 
iLutz et all (|2010h and IShao et~afl (|2010l ). Their steeper 
power law, with a slope of 0.8, for luminous quasars is 



not directly comparable with our result as it does not 
separate out the independent effect of redshift on Lir. 
However, we would also note that the quasar part of 
this re lationship is based on samples presented bv lNetzerl 
(|2009l . and earlier works) , which ignore AGN without de- 
tections in the far-infrared or mm, so that the observed 
correlation may be partly due to far-infrared flux limits. 

To investigate whether our fit could be affected by a 
redshift-dependence in 8, we perform the maximum like- 
lihood fitting again, using objects in restricted redshift 
ranges where we have enough objects to make a robust 
measurement, and marginalise over all parameters except 
0. The results of this are shown in Figure [5] the derived 
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values of 9 are not significantly different from each other 
for redshifts 0.5 < z < 2.5, although the small number 
of objects in the redshift bins means we cannot make a 
strong statement for or against evolution of 9 as a func- 
tion of redshift. However, the 0.5 < z < 1 range shows a 
detection of non-zero 9 at better th an the 2-a level, which 
is com plementary to the result of Hat ziminaoglou et al.l 
(2010). While we do not have enough data to detect any 
significant decline in 9 as a function of redshift, the data- 
po ints in Figure [5] arc consistent with the relation found 
bv lSerieant et aljlfeoich . and with our smaller H-ATLAS- 
only dataset we cannot rule out the existence of such a 
weak redshift dependence. 

Since the far-infrared lu minosity Ltr has been show n 
to be due to star formation |Hatziminaoglou et al.ll20l"oT ). 
the physical interpretation of these results is that star- 



formation rate is higher both as redshift increases and 
as qu asar luminosity increases. As was also suggested by 
SH09J, and as has been incorporated into the "quasar- 
mode" feedback r ecipes of semi-analytic models (e.g. 
ICroton et afll2006h . one of the simplest ways to inter- 
pret this is if both black hole accretion and star for- 
mation depend on a common supply of cold g 
might be suppli ed during a galaxy-galaxy merger or cold- 
accretion event (Dckcl et al. 2009; Dijkstra & Loeb 2009; 



ISmith fc Jarvisl l2007t ISmith et al.l 120081 1. The implica- 
tion of the non-unity-slope is that the dependence of ei- 
ther star-formation or accretion luminosity (or both) on 
the gas supply cannot simply be linear with gas mass, as 
one might naively suppose, but is instead more complex. 

For both star formation and black-hole accretion, 
such gas must be sufficiently cool that pressure support 
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1.0 r 



0.5 



-0.5 



during mergers or cold accretion, or heating of gas via 
feedback processes. 

With the full 550 square degrees of the H-ATLAS 
survey, we will be able to substantially increase the sam- 
ple of quasars and improve these constraints. The larger 
area will allow us to obtain representative samples of the 
objects with the lowest density on the sky, i.e. quasars 
with high luminosities at low redshift, which will increase 
the sensitivity of our test at low redshifts and permit us to 
remove the assumption of power-law correlations and test 
more general models. More importantly, a larger sample 
of objects at all redshifts will make it possible to measure 
the temperatures of dust in quasar hosts as a function of 
Lqso and redshift by stacking the far-infrared fluxes; this 
in turn will remove the largest source of uncertainty in 
the relationship between Lib and z. 



-1.0 1 i i i i i i i i i i i i i i i i i i i I 

0.5 1.0 1.5 2.0 2.5 

Redshift 

Figure 5. Best-fitting values of based on maximum- 
likelihood fits using data in the redshift ranges denoted by 
the widths of the errorbars. Black, solid errorbars are 68.3 per 
cent confidence limits; blue, dashed errorbars are 95.4 per cent 
confidence limits. 



does not prevent it from being accreted by the AGN or 
from collapsing under its own gravity to form stars, so an 
alternative way to make this connection is via a mecha- 
nism (such as AGN or supernova feedback) which is ca- 
pable of heating gas at all scales within a galaxy. 



5 CONCLUSIONS 

Using far-infrared measurements (at 250, 350 and 
500 urn) of SDSS and 2SLAQ quasars in the H-ATLAS 
science demonstration field, we find evidence for corre- 
lations between the infrared luminosities of quasar hosts 
with both redshift and quasar luminosity. While the best- 
fit power-law index for the correlation with redshift can 
vary significantly depending on the assumed tempera- 
ture of the dust in the quasar hosts, the correlation with 
quasar luminosity is robust to this effect, and is non- 
zero at the 2.5-cr level when we combine all objects in a 
maximum- likelihood fit. 

We find a power-law index 6 = 0.22 ± 0.08 for the 
relationship between Lib and L qso . Our data show that 
9 is non-zero at the 2-a level in the 0.5 < z < 1 redshift 
range alone, but we have too few objects to find evidence 
for evolution in 8 over the range 0.5 < z < 2.5. 

We find evidence for a large intrinsic dispersion in the 
redshift dependence, but no evidence for intrinsic disper- 
sion in the correlation between Lqso and Lm, suggesting 
that the latter is due to a direct physical connection be- 
tween star formation and black hole accretion. One pos- 
sible interpretation of this is that both the quasar activity 
and star formation are dependent on the same reservoir 
of cold gas, and are thus both affected by influx of gas 
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